Molecular Dynamics Simulation of Semiflexible Polyampholyte Brushes - The Effect 

of Charged Monomers Sequence 
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Planar brushes formed by end-grafted semiflexible polyampholyte chains, each chain containing 
equal number of positively and negatively charged monomers is studied using molecular dynamics 
simulations. Keeping the length of the chains fixed, dependence of the average brush thickness and 
equilibrium statistics of the brush conformations on the grafting density and the salt concentration 
are obtained with various sequences of charged monomers. When similarly charged monomers of the 
chains are arranged in longer blocks, the average brush thickness is smaller and dependence of brush 
properties on the grafting density and the salt concentration is stronger. With such long blocks of 
similarly charged monomers, the anchored chains bond to each other in the vicinity of the grafting 
surface at low grafting densities and buckle toward the grafting surface at high grafting densities. 

PACS numbers: 



I. INTRODUCTION 

Polymers carrying ionizable groups dissolved in a po- 
lar solvent dissociate into charged polymer chains and 
counter- ions (ions of opposite charge). Depending on 
acidic or basic property of ionizable monomers, charged 
polymers can be classified into polyelectrolytes contain- 
ing a single sign of charged monomers and polyam- 
pholytes bearing both positive and negative charges. 
These macromolecules are often water-soluble and have 
numerous applications in industry and medicine. Many 
biological macromolecules such as DNA, RNA, and pro- 
teins are charged polymers. 

Charged monomers of different sign can be distributed 
randomly along a polyampholyte chain or charges of one 
sign can be arranged in long blocks. With excess of one 
type of charges, polyampholyte solutions exhibit proper- 
ties similar to those of polyelectrolyte solutions. With 
the same ratio of positive and negative charges on the 
chains, the solution behavior depends noticeably on the 
charge sequence. For example, it has been shown that 
the sequence of charged amino acids (charge distribu- 
tion) along ionically complementary peptides affect the 
aggregation behavior and self-assembling process in the 
solution of such peptides [U, Also, using Monte Carlo 
simulations it has been shown that charged monomers 
sequence of neutral block-polyampholytes affect their ad- 
sorption properties to a charged surface 

The present understanding of, in particular, random 
polyampholytes in solution and their interaction with 
surfaces and polyelectrolytes has recently been reviewed 
In addition, theoretical and computer simulation 
studies of single diblock polyampholytes d, @, and 
diblock polyampholytes in solution [1, Q , have been per- 
formed. The value of salt concentration in the solu- 
tion is an important parameter for tuning the structural 



properties of polyampholyte chains. With increasing the 
salt concentration, the coil size and the viscosity in- 
crease in the solution of nearly charge-balanced polyam- 
pholyte chains and decrease in the case of the solution 
of polyampholytes with a larg e net charge (the polyelec- 
trolyte regime) 0, [H 113, M- 

The properties of the system of polymers anchored on 
a surface are of great interest both in industrial and bio- 
logical applications and academic research. When there 
is sufficiently strong repulsion between the polymers, the 
chains become stretched and the structure obtained is 
known as a polymer brush. Brushes at planar and curved 
surfaces formed by grafted homopolymers have exten- 
sively been investigated by various theoretical methods 

In the case of polyelectrolytes grafted onto a surface 
(polyelectrolyte brush), the repulsion of electrostatic ori- 
gin between the chains can be considerable even at low 
grafting densities, making it easy the system to access the 
brush regime. The use of charged polymers instead of 
uncharged polymers introduces additional length scales 
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FIG. 1: (Color online) Schematic structures of four different 
chains containing different sequences of charged monomers. 
Neutral, positively and negatively charged monomers are 
shown by white, red (dark) and blue (black) spheres re- 
spectively. The chain containing only negatively charged 
monomers is a polyelectrolyte chain, PE, and three other 
chains PA 11, PA33 and PA66 are polyampholytes consisting 
of alternating blocks of similarly charged monomers with 1, 3 
and 6 monomers in each block respectively. 
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FIG. 2: (Color online) Typical equilibrium configurations of 
brushes formed by a)PE, b)PAll, c)PA33, and d)PA66 chains 
at grafting density paO"^ = 0.1. In the upper plans, counte- 
rions are not shown and periodic boundary condition is re- 
moved for clarity. As it can be seen, in the brushes of PA33 
and PA66, the chains are buckled toward the grafting surface. 



such as Bjerrum length and Debye screening length to 
the system. Polyelectrolyte brushes also have been sub- 
jected to extensive investigations involving both theoret- 
ical [H, M, 

lation methods 




and computer simu- 
2i|. At high enough 
grafting densities and charge fractions of polyelectrolyte 
chains, most of counterions are trapped inside the poly- 
electrolyte brush and competition between osmotic pres- 
sure of the counterions and elasticity of the chains de- 
termines the brush thickness. This regime of a polyelec- 
trolyte brush is known as the osmotic regime in which 
some theoretical scaling methods predict no dependence 
of the brush thickness to the grafting density [11, H^. 
However, other scaling method that takes into account 
the excluded volume effects and nonlinear elasticity of 
polyelectrolyte chains predicts a linear dependence of the 
brush thickness on the grafting densit y a, nd is in agree- 
ment with experiment and simulation |23l. [2^ [25j. Also, 
it has been shown that diffusion of a fraction of coun- 
terions outside the polyelectrolyte brush leads to a log- 
arithmic dependence of the average brush thickness on 
the grafting density [2^ . 

Brushes formed by grafted diblock polyampholytcs 
have also been investigated by lattice mean field modeling 
[3^, 'sT'l and computer simulation [33| . The effect of chain 
stiffness, charge density, and grafting density on spheri- 
cal brushes of diblock polyampholytes and interaction be- 
tween colloids with grafted diblock polyampholytes have 
been studied using Monte Carlo simulations [33, ,34] . De- 
spite the case of an osmotic polyelectrolyte brush that 
counterions are trapped inside the brush and have a cru- 
cial rule in the equilibrium brush thickness, in the case of 
a polyampholyte brush most of counterions are outside 
and equilibrium brush thickness is mainly determined by 



the chains properties. 

A brush of semiflexible charged polymers is a dense as- 
sembly of such polymers in which the interplay between 
electrostatic correlations, entropic effects and the bend- 
ing elasticity of the chains can bring about a variety of 
equilibrium properties. The strength of the electrostatic 
correlations in a polyampholyte brush depends on the se- 
quence of charged monomers along the chains. With sim- 
ilarly charged monomers being arranged in long blocks, 
the electrostatic correlations are of great importance in 
the system. It seems that for a brush formed by semiflex- 
ible polyampholyte chains, each chain containing equal 
number of positive and negative monomers, studying the 
effect of monomers sequence on equilibrium properties of 
the brush could be an step toward better understanding 
of such system s. 

In this paper, we study planar brushes of semiflexi- 
ble polyampholytes using extensive molecular dynamics 
(MD) simulations in a wide range of the grafting density 
and at various salt concentrations. In our simulations, 
each end-grafted polyampholyte chain contains equal 
number of positively and negatively charged monomers 
(isoelectric condition). We investigate the effect of the 
sequence of charged monomers along the chains on equi- 
librium properties of the brush. Charged monomers of 
each sign are arranged in blocks of equal length and 
monomers of opposite sign form regular alternation of 
positive and negative blocks (see Fig. [T] b, c and d). 
Different charge sequences in our simulations correspond 
to different lengths of above mentioned blocks. We also 
consider the case of polyelectrolyte brush in our studies 
in which all the charged monomers of the chains are of 
the same sign for comparison. 

We find that at each value of the grafting density, the 
brush formed by polyampholytes with longer blocks of 
similarly charged monomers has smaller average thick- 
ness, i.e., the lowest value of the average thickness cor- 
responds to the brush of diblock polyampholytes. The 
values of the grafting density that we use in our simula- 
tions correspond to the osmotic regime of the polyelec- 
trolyte brush. In this regime, most of the counterions are 
contained inside the polyelectrolyte brush and the aver- 
age brush thickness has a weak linear dependence on the 
grafting density. Dependence of the average thickness 
of polyelectrolyte brush that we obtain from our simula- 
tions is in agreement with nonlinear theory of ref. [23| . 
In the case of the polyampholyte brush, we find that 
the average thickness considerably varies as a function of 
the grafting density and with longer blocks of similarly 
charged monomers, this dependence is stronger. Coun- 
terions of a polyampholyte brush are not trapped inside 
the brush and dependence of the average thickness on 
the grafting density originates from a different mecha- 
nism than that of a polyelectrolyte brush (see the text). 
We also study the equilibrium statistics of conformations 
of the chains at various grafting densities and salt con- 
centrations both for polyampholyte and polyelectrolyte 
brushes. We find that in the case of polyampholyte 
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FIG. 3: (Color online) Average thickness of brushes formed 
by PE, PAll, PA33 and PA66 chains versus dimensionless 
grafting density paO^ with no added salt. Dotted line is the 
prediction of the scaling theory of ref. [2^ with o^jf ~ 1.02(t 
and dashed line is a linear fit to our data for the brush of 
PA66 chains. As it can be seen dependence of the average 
thickness of polyampholyte brushes with longer blocks of sim- 
ilarly charged monomers on the grafting density is stronger. 
The size of the symbols corresponds to the size of error bars. 



brushes with long blocks of similarly charged monomers, 
at low grafting densities the inter-chain electrostatic in- 
teractions link the chains to each other in the vicinity of 
the grafting surface and cause the brush to collapse. At 
high grafting densities, although the excluded volume in- 
teractions tend to increase the brush thickness, the elec- 
trostatic interaction dominate over the bending elastic- 
ity of the chains and lead them to buckle and the brush 
thickness to decrease. Also it is observed that with long 
blocks of similarly charged monomers along the chains, 
the average thickness of the polyampholyte brush is an 
increasing function of the salt concentration. 

The rest of the paper is organized as follows. In Sec. 
we describe our model and details of molecular dynamics 
simulations. We present results of the simulations in de- 
tail and concluding remarks in Sec. IIIII Finally, in Sec. 
|V]we summarize the paper. 



II. THE MODEL AND THE SIMULATION 
METHOD 

In our simulations which areperformed with the MD 
simulation package ESPResSo [3a| , each brush is modeled 
by M semiflexible polyelectrolyte/polyampholyte bead- 
spring chains of length N (N spherical monomers) which 
are end-grafted onto an uncharged surface at z — 0. 
Short-range repulsive interaction which is described by 
a shifted Lennard- Jones potential. 



\} if r < Tc, 
if r > Tr, 



(1) 



is considered between monomers in which e and a are 
the usual Lennard-Jones parameters and the cutoff ra- 
dius is Tc — 2^/^(7. Neighboring beads along each chain 
are bonded to each other by a FENE (finite extensible 
nonlinear elastic) potential |36| , 



Ubond{r) 



ifr>i?o, 



(2) 



with bond strength fctond = SOe/ct^ and maximum bond 
length i?o = 1.5(7. The chains are semiflexible and their 
elasticity is modeled by a bond angle potential, 



Ube7id{r) = fche„<i(l - COS 6*), 



(3) 



in which 9 is the angle between two successive bond vec- 
tors of a chain. In our simulations kbend = ^bksT which 



means that Ip = 25a 



Lc in which Ip and Lc are the 



persistence length and the contour length of the chains 
respectively. The simulation box is of volume L x L x Lz 
in which L is the box width in x and y directions and 
Lz is its height in z direction and the grafting density 
is given by pa — M / L? . The positions of anchored 
monomers which are fixed during the simulation, form 

— 1/2 

an square lattice with lattice spacing d — pa on the 
grafting surface {x — y plane). All particles interact re- 
pulsively with the grafting surface at short distances with 
the shifted Lennard-Jones potential introduced in Eq. [1] 
In addition, a similar repulsive potential is applied at the 
top boundary of the simulation box and in our simula- 
tions Lz — 2Na. We model positive and negative ions of 
monovalent salt by equal number of spherical Lennard- 
Jones particles of diameter a with charges e and — e re- 
spectively. All the charged particles interact with the 
Coulomb interaction 



uc{r) = ksTqigj- 



(4) 



in which Qi and Qj are charges of particles i and j in 
units of elementary charge e and r is separation between 
them. The Bjerrum length, Ib, which determines the 
strength of the Coulomb interaction relative to the ther- 
mal energy, fc^T, is given by Ib = e^/ekBT, where e is 
the dielectric constant of the solvent. Ib — "ia in our sim- 
ulations. We apply periodic boundary conditions only in 
two dimensions (x and y). To calculate Coulomb forces 
and energies, we use the so-called MMM technique in- 
troduced by Strebel and Sperb [37l| and modified for lat- 
erally periodic systems {MMM2D) by Arnold and Holm 
[sst . The temperature in our simulations is kept fixed at 
kBT = 1.2e using a Langevin thermostat. 

In our simulations we consider brushes containing M ~ 
25 anchored chains, each chain consisting of iV = 24 
monomers, fN of them are charged and / = 5- We 
also consider M^E. counterions to neutralize the chains 
charge. The following sequences of charged monomers 
along the chains are considered to model polyelectrolyte 
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Csa^ = 0.22. In the beginning of each simulation, all 
of the chains are straight and perpendicular to the graft- 
ing surface and all the ions are randomly distributed in- 
side the simulation box. We equilibrate the system for 
1.6 X 10^ MD time steps which is enough for all values of 
the grafting density mentioned above and then calculate 
thermal averages over 1500 independent configurations of 
the system selected from 2.25 x 10® additional MD steps 
after equilibration. MD time step in our simulations is 



T = O.OIto in which tq = y is the MD time scale 
and m is the mass of the particles. 

We calculate the average brush thickness which can 
be measured by taking the first moment of the monomer 
density profile 



/q°° Zprn{z)dz 
Pm{z)dz ' 



(5) 



in which pmiz) is the number density of monomers as 
a function of the distance from the grafting surface. 
For better understanding of the statistics of the chains 
conformations, we calculate the histogram of the mean 
end-to-end distance of the chains, P{R), in which B = 

jjJ2iLi l-Ril and Ri is the end-to-end vector of chain i. 
We also calculate the histogram of the average distance 
of the end monomers of the chains from the grafting sur- 
face, P{zend), in which Zend = jjJ^iLi^i ^ud Zi is the 
z component of the end monomer of chain i. Also, as a 
measure of the lateral fluctuations of the chains we define 



hat as hat 



in which 



II 



- Ri.z\ is 



the magnitude of the lateral component of Ri. 



ST 



rh 



FIG. 4: The histogram of the average end-to-end distance of 
the chains for brushes formed by PE, PA 11, PA33 and PA66 
chains at grafting densities a) paO"^ = 0.1, h)paU^ = 0.06 and 



cj pair 



0.02. 



and three different polyampholyte brushes respectively: 

PE (. .. ...) PAll {...+ - + - + -...), PA33 

(...+ + + ... )and PA66 (...+ + + + + + -- 

- - - - ...). Each chain consist of an alternating se- 
quence of charged and neutral monomers and the neutral 
monomers are not shown for simplicity (see Fig. [T]). We 
do simulations of the brushes formed by four above men- 
tioned different chains at dimensionless grafting densities 
p^cr^ = 0.02,0.04,0.06,0.08,0.10 and at various mono- 
valent salt concentrations changing from c^a^ = to 



III. RESULTS AND CONCLUDING REMARKS 

We calculate the average thickness of the brushes 
formed by polyelectrolyte chain, PE, and three different 
polyampholyte chains PAll, PA33 and PA66 at different 
grafting densities. In Fig. [3] the average brush thick- 
ness versus dimensionless grafting density pa<J^ with no 
added salt is shown. As it can be seen in this figure, the 
average thickness of the polyelectrolyte brush has a week 
dependence on the grafting density. With grafting densi- 
ties that we use in our simulations, the Gouy-Chapman 
length of the polyelectrolyte brush which is defined as 
Xg = {"^T^PafNls)^^ , is quite smaller than the length 
of the chains {Xg < cr) and most of the counterions are 
contained inside the brush (see Fig. [T]a). In this regime 
of polyelectrolyte brushes, the osmotic pressure of the 
counterions determines the thickness of the brush and the 
week dependence of the brush thickness on the grafting 
density is expected. In fig. [3] the prediction of the scaling 
theory of ref. [I^] (see Eq. (28) of [23]) corresponding 
to our parameters is shown for comparison. As it can be 
seen, this scaling theory describes well the dependence 
of the average thickness of semiflexible osmotic polyelec- 
trolyte brush on the grafting density obtained from our 
simulations. 
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FIG. 5: The histogram of the average distance of the end 
monomers of the chains from the grafting surface for brushes 
formed by PE, PAll, PA33 and PA66 chains at grafting den- 
sities a) paCT^ = 0.1, b)paO"^ = 0.06 and c) paa^ = 0.02. 



In the case of the brushes of polyampholyte chains, the 
net charge of the brush layer is zero and most of the coun- 
terions are outside the brush. The average thickness in 
this case has a considerable dependence on the grafting 
density and the interaction between the anchored chains. 
At each value of the grafting density, brush thickness de- 
creases with increasing the length of the blocks of simi- 
larly charged monomers. Also, as it can be seen in Fig. O 
the average thickness of the brushes formed by the chains 
with longer blocks of similarly charged monomers has 



stronger dependence on the grafting density (see grow- 
ing slope of -^^i^ versus paC^ from PAll to PA66). 

To describe such dependence of the brushes thickness 
on the grafting density it is instructive to study the equi- 
librium statistics of conformations of the chains. Ac- 
cordingly, we monitor the histograms P{R) and P{zend) 
which are defined in Sec. |TI] in our simulations. In Figs. 
mandOthe histograms P{R) and P{zend) at three differ- 
ent grafting densities are shown for brushes of polyelec- 
trolyte and three different polyampholyte chains. It can 
be understood from Figs. 2] and [5] that decreasing the av- 
erage thickness of a polyampholyte brush with increasing 
the length of the blocks of similarly charged monomers 
originates from two different mechanisms at low and high 
grafting densities. At low grafting densities, the anchored 
chains both in polyelectrolyte brush and in three differ- 
ent polyampholyte brushes have stretched conformations 
and the average end-to-end distance, R, has its maxi- 
mum value in most of equilibrium configurations (see 
Fig. [4]c). In Fig. [5] c, the histogram P{zend) shows 
that at such low values of the grafting density the poly- 
electrolyte brush is aligned and the PE chains are mostly 
perpendicular to the grafting surface. In this figure, it 
can be seen that in the case of polyampholyte brushes 
with long blocks of similarly charged monomers, 2 = 
part of P{zend) IS nonzero and increases with increas- 
ing the length of blocks of similarly charged monomers. 
From these two histograms one can conclude that in the 
brushes of polyampholytes with long blocks of similarly 
charged monomers, PA33 and PA66, the chains have rod- 
like conformations and fluctuate in the vicinity of the 
grafting surface. Snapshots of the brushes of these chains 
(not shown here) show that the chains are linked to each 
other as doublets near the grafting surface because of the 
strong electrostatic attractions. In the case of the brush 
of PAll chains, the end monomers of the chains have a 
uniform distribution inside the brush. In this case, elec- 
trostatic neutrality is satisfied in smaller length scales 
and the electrostatic correlations between the chains are 
very weak. 

At high grafting densities, the histogram of the aver- 
age end-to-end distance of the chains shows that in the 
brushes of polyampholyte chains consisting of long blocks 
of similarly charged monomers, despite the case of low 
grafting density, the chains tend to have buckled con- 
formations. The histogram of the average z component 
of the end monomers, P(ze„d), shows two maximums 
around z = and z = L in the case of the brush of 
diblock PA66 chains and three maximums in the case of 
the brush of PA33 chains. The energy needed for buck- 
ling a semiflexible chain (Lc ~ Ip), Ef,, is of the order 
of thermal energy, ksT. In dense brushes of the chains 
PA33 and PA66, blocks of similarly charged monomers of 
the chains form highly charged layers of monomers which 
produce strong electrostatic field. When the electrostatic 
energy of similarly charged monomers of the end blocks 
of the chains due to the electric field of the layer of op- 
posite charge exceeds the bending energy, Eb ^ ksT^ 
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FIG. 6: Salt concentration dependence of the average thick- 
ness of brushes of PAll and PA66 chains at grafting density 
PaO"^ — 0.08. As it can be seen, the thickness of the brush 
formed by the chains with long blocks of similarly charged 
monomers is an increasing function of the salt concentration. 
The size of the symbols corresponds to the size of error bars. 



the chains buckle toward the grafting surface. For ex- 
ample, in the case of the brush of diblock PA66 chains, 
the electrostatic energy of positively charged monomers 
of each chain roughly can be written as E^i — ^-^E^ 
in which E — Airpae/e = AirkBTlBPa/^ is approximately 
the electric field of the layer of anchored blocks of neg- 
atively charged monomers. With parameters that we 
use in our simulations, at high grafting densities of the 
brushes formed by PA33 and PA66 chains, Eei > ksT 
and buckling of the chains is expected. In the case of the 
brush of PA66 chains, two maximums in the histogram 
P{zend) shows that buckled chains and those that are 
straight and perpendicular to the grafting surface coex- 
ist. In the brush of PA33 chains, there are two layers of 
negatively charged monomers and considering the chains 
of straight configuration, three maximums seen in the 
histogram P{zend) is expected. 

The effect of added salt on the brushes of polyam- 
pholyte chains is also studied using explicit monovalent 
salt ions. We observe that properties of the brushes 
formed by the chains containing longer blocks of simi- 
larly charged monomers are more sensitive to the salt 
concentration. The average thickness of these brushes at 
all grafting densities increase with increasing the salt con- 
centration and reach to its maximum value (see Fig. 
This behavior is completely different from that of a poly- 
electrolyte brush which its thickness is known that de- 
creases with increasing the salt concentration. In the case 
of a polyelectrolyte brush the electrostatic repulsion be- 
tween the chains tend to swell the brush and the screen- 
ing effect of added salt decreases the average brush thick- 
ness. However, in the case of the brush of polyampholyte 
chains with long blocks of similarly charged monomers, 
the electrostatic interactions between the chains are at- 
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FIG. 7: Histograms a)P(2end) and h)P(R) for the brush 
formed by PA66 chains at different salt concentrations. The 
grafting density is paO^ = 0.08. 



tractive and tend to decrease the brush thickness. In 
these brushes, electrostatic screening of the salt ions in- 
creases the average thickness. It also can be seen in Fig. 
[6] that the added salt doesn't change the average thick- 
ness of the brush of PAll chains in which because of local 
electrostatic neutrality the electrostatic correlations are 
negligible. In Fig. [7] the equilibrium histograms P{R) 
and P{zend) for the brush of PA66 chains are shown at 
different salt concentrations. It can be seen that with 
adding more salt to the system, the chains become more 
straight (see part b of the figure) and z = part of the 
histogram P{zend) vanishes which means that the brush 
becomes more aligned. 

To estimate the importance of finite size effect in our 
simulations, we calculated equilibrium average of liat (de- 
fined in Sec. [n| for the brush of PA66 chains at highest 
grafting density PaC^ — 0.1 and found that {Uat) — 11. Tcr. 
The fact that the value of {hat) is less than the lateral box 
length, L ~ 16cr, means that the chains doesn't overlap 
with their own images in our simulations. We also re- 
peated the simulation of the PA66 brush at grafting den- 
sity pa<J^ = 0.1 with M = 64 chains (instead of M = 25) 
and didn't obtain different results. 
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IV. DISCUSSION 

Osmotic pressure of trapped counterions inside poly- 
electrolyte brush in the osmotic regime is one of the 
main factors determining the equilibrium brush thick- 
ness. Nonetheless, because of the electrostatic neutral- 
ity of the chains, Gouy-Chapman length diverges in the 
case of polyampholyte brushes considered in this paper 
and most of counterions are outside the brush layer (see 
Fig. [2]b, c and d). These counterions have no rule in the 
equilibrium brush thickness. In this case the brush thick- 
ness is resulted from the interplay between inter- and 
intra-chain electrostatic attraction of oppositely charged 
blocks, inter-chain excluded volume interactions and the 
bending elasticity of the chains. Hence, dependence of 
the average thickness of such polyampholyte brushes on 
the grafting density is of quite different mechanism. How- 
ever, it is interesting that this dependence is still linear 
similar to that of a polyelectrolyte brush (see Fig. [3]). 
This dependence in brushes of polyampholyte chains with 
long blocks of similarly charged monomers is appreciably 
stronger than the case of polyelectrolyte brush (for ex- 
ample in Fig. O the slopes of dotted and dashed lines 
are 16.7 and 40.2 respectively). An important parameter 
in a brush of semiflexible polyampholytes is the bending 
rigidity of the chains, kbend- With decreasing the value of 
kbend^ electrostatic attractions dominate over the chains 
stiffness and in competition with excluded volume effects 
determine the average brush thickness. The opposite ex- 
treme is the case of the brush of rod-like polyampholytes 
in which electrostatic correlations can not buckle the 
chains. So, different regimes corresponding to different 



values of ktend should be investigated (39j . 

V. SUMMARY 

In summary, we have used molecular dynamics simu- 
lations to study planar brushes formed by semiflexible 
polyelectrolytes and polyampholyte chains with differ- 
ent sequences of charged monomers at various grafting 
densities and salt concentrations. It has been shown 
that at grafting densities corresponding to the osmotic 
regime of the polyelectrolyte brush, the average brush 
thickness is a weak linear function of the grafting density 
in agreement with predictions of ref. [23|. The aver- 
age thickness and equilibrium properties of the polyam- 
pholyte brushes have considerable dependence on the 
grafting density and the salt concentration. This depen- 
dence is stronger for brushes of polyampholyte chains 
containing longer blocks of similarly charged monomers. 
In brushes of polyampholyte chains with long blocks of 
similarly charged monomers, at low grafting densities the 
electrostatic attractions link the chains to each other 
in the vicinity of the grafting surface and collapse the 
brush. At high grafting densities, the electrostatic corre- 
lations dominate over the bending rigidity of the chains 
and cause them to buckle toward the grafting surface. 
Despite the case of polyelectrolyte brushes, the average 
thickness of polyampholyte brushes with long blocks of 
similarly charged monomers increases with increasing the 
salt concentration. 

We are grateful to A. Naji, R. Golestanian and F. 
Mohammad- Rafiec for useful comments and discussions. 
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